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We report measurements of charmless hadronic decays into the Tv^n~ K'^-k~ final state. The 
analysis uses a sample of 657 x 10® BB pairs collected with the Belle detector at the KEKB 
asymmetric-energy e'^e~ collider at the T(4S) resonance. The decay B" — > p^K^-K~ is observed 
for the first time; the significance is 5.0o" and the corresponding partial branching fraction for 
Mifvr G (0.75, 1.20) GeV/c^ is [2.8 ± 0.5(stat) ± 0.5(syst)] x lO"**. We also obtain the first evidence 
for _B° — > /o(980)A"''7r~ with 3.5(t significance and for — > ■k'^-k~K*^ with 4.5cr significance. 
For the two-body decays p^K*^ and _B° /o(980)A'*'', the significances are 2.7(t and 2.5cr, 

respectively, and the upper limits on the branching fractions are 3.4 x 10"*^ and 2.2 x 10~® at 90% 
confidence level. 



PACS numbers: 11.30.Er, 13.25.Hw, 14.40.Nd 



In the Standard Model (SM), charmless hadronic B 
meson decays occur mainly via two processes: (i) b 
sqq transitions mediated by penguin diagrams, and (ii) 
b uW* transitions mediated by tree diagrams. These 
diagrams as they pertain to B^ p'^K*'^ for ex- 
ample, are shown in Fig. [TJ Both of these processes are 
suppressed relative to the more common b — > cW* decays 
due to either (i) the one-loop structure or (ii) the small 
ratio of CKM matrix elements iKtb/Vcbl, respectively. 
Because of this suppression, these decays are especially 
sensitive to non-SM contributions 

There have been several puzzling results from mea- 
surements of charmless hadronic B decays. For exam- 
ple, B decays to K'^tt' and K^tt^ show different pat- 
terns of direct CP violation Q, which are inconsistent 
with naive SM expectations. It has been suggested Q 
that vector-vector (VV) final states with the same quark 
combinations, e.g. B — > pK* may give insights to this 
puzzle, as any difference between Kt: and their VV coun- 
terparts will be mainly hadronic. In addition, charmless 
B decays to VV final states show intriguing results in 




the final-state polarizations. The decays B <f)K* and 
B pK* , both occurring mostly via the b ^ s pen- 
guin process, are found to have large transverse polariza- 
tions @, 0, Si , in contrast to the expectation from factor- 
ization. On the other hand, iJ+f"' which is 
mostly a 6 — > uW* tree-diagram process, is almost fully 
polarized longitudinally There have been theoreti- 
cal studies of these modes, in part focusing on the 
final-state polarizations within and beyond the SM. 

One difficulty in measuring charmless B VV de- 
cays, however, is that non-resonant decays to the same 
final state can be a significant background Such 



FIG. 1: Feynman diagrams for charmless hadronic B decays 
pertaining to B" — > p'^K*'^: (a) 6 ^ s penguin diagram, (b) 
b ^ u tree diagram. 



non-resonant decays may have different decay properties, 
e.g., different polarization of the vector mesons. While 
there arc several experimental studies of B^ decays to 
the 7r+7r~X+7r~ final state 0, there is no experi- 
mental information on the non-resonant components of 
these final states. 

In this paper, wc analyze charmless hadronic de- 
cays of B^ to the Tr+TT^if+Tr^ final state. We search 
for two-body final states such as — > p^K*^ and 
/o(980)if*°, and also for three-body states p°K+Tr^, 
/o(980)is:+7r- and tt+tt-R*'^, where the tt+tt- or 

K'^TT^ pairs are non-resonant. A comprehensive under- 
standing of these decays with a clear distinction between 
non-resonant and two-body resonant decays would ad- 
vance our understanding of strong and weak interaction 
dynamics. 

We use a data sample containing 657 x 10^ BB 
pair events collected with the Belle detector [iJl at 
the KEKB [l^ asymmetric-energy e+e" collider (3.5 
on 8 GeV), operating at the T(4S') resonance. To re- 
construct B'^ — > 7r"'"7r~_ft'"'"7r~ decays including the in- 
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termediate states — > tt+tt", /o(980) 'k'^'k~ , and 
_^ K~^Tr~ , we select four charged tracks of which two 
are positively charged and two are negatively charged. 
Each track is required to originate within 5.0 cm of the in- 
teraction point (IP) along the beam direction, and within 
0.2 cm of the IP in the transverse plane 01 ■ We also re- 
quire that the transverse momentum of each track be 
larger than 0.1 GeV/c Tracks identified as electrons 
are rejected. We identify charged kaons and pions by 
combining particle identification (PID) information ob- 
tained from the central drift chamber, the time-of-flight 
system, and the aerogel Cherenkov counters fl^ . 

Signal candidates are selected for further analysis 
based on four kinematic variables: the 7r"''7r~ and K~^'!t~ 
invariant masses (Mtttt and Mk-k), the energy difference 
AE = Eb — -Ebcam, and the beam-energy-constraincd 
mass Mbc = V^bcam ~ Pb^ where -Eboam is the beam 
energy and Eb and pb are the energy and momen- 
tum, respectively, of the candidate i?-meson. These vari- 
ables are all evaluated in the T(4S') center-of-mass (CM) 
frame. We retain events satisfying \AE\ < 0.1 GeV, 
5.24 GcV/c2 < Mbc < 5.29 GeV/c^, 0.55 GeV/c^ < 
M^^ < 1.20 GcV/c2 and 0.75 GeV/c^ < Mk^ < 
1.20 GcV/c^. To optimize the background suppression 
criteria, tighter "signal regions" are defined for Mbc and 
AE: 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and \AE\ < 
0.045 GeV. The fraction of events having multiple can- 
didates is approximately 20%. For multiple-candidate 
events, we select the candidate decay having the smallest 
from the B vertex fit. Given a set of four particles, 
TT^TT^ K^Tr~ , two combinations of (tt+tt") and (K~^tt~) 
may lie inside the selected mass ranges. In this case, 
which occurs in less than 1% of signal decays, we pair 
the higher-momentum 7r~ with the n'^ . 

The dominant source of background is from continuum 



qq events (g 



d, s, c). These events are dis- 



tinguished from the signal by their event shape. Since B 
mesons are spinless and produced nearly at rest in the 
GM frame, their daughter particles are distributed al- 
most isotropically. On the other hand, continuum events 
usually produce two back-to-back jets in the CM frame. 
We use Monte Carlo (MC) simulated [l3| signal events 
and sideband data (5.20 GcV/c^ < Mbc < 5.26 GcV/c^) 
for optimizing the continuum suppression requirements. 
First wc form a Fisher discriminant J- based on a set 
of modified Fox- Wolfram moments [T^j. These moments 
are uncorrclated with the four kinematical variables men- 
tioned above. Two more variables are used for continuum 
suppression: cos 9b, the cosine of the polar angle of the 
B flight direction in the CM frame; and Az, the displace- 
ment along the beam direction between the vertex of the 
signal B and that of the other B in the event. Likeli- 
hood functions for signal (Cs) and continuum background 
{Cqq) are formed from products of the probability den- 
sity functions (PDFs) for T, cos9b, and Az. These arc 



To obtain improved continuum suppression, we optimize 
the requirement on TZqq as a function of flavor-tagging 
information from the accompanying B meson. The Belle 
flavor-tagging algorithm [l^ yields the 6-flavor variable q 
(=±1), and the quality variable r. The latter ranges from 
zero for no flavor discrimination to one for unambiguous 
flavor assignment. We optimize the TZqq requirement in- 
dependently in six bins of qr. For B'^ — s- p^K*^, for 
example, the optimized Ttqq requirements remove 99% of 
the qq background while retaining 42% of the signal. 

B decays to a charm meson (_D" or £)'*)+) and multiple 
pions constitute a signiflcant background that exhibits 
peaking behavior in Mbc and AE similar to that of the 
signal. To eliminate this background, wc veto candidates 
that have a Kmr, Kit, or tttt invariant mass consistent 
with a K-TT+TT+, £>" K-7r+, and D" 

7r~7r+ decay, respectively. 

The signal yields are obtained from a four-dimensional 
extended unbinned maximum-likelihood [lOl fit (4D fit) 
to Mbc, AE, AfTTTT and Mktt- The likelihood function is 



(1) 



1=1 j 



combined into a likelihood ratio TZa 



where Yj is the yield of the j-th component, V^j is 
the PDF value for the j-th component of the i-th 
event, and i runs over all events in the fit region {N). 
We include 13 components in Eq. (P): B° decays to 
p°K*°, /o(980)if*^ and f2{l21Q)K*^; the non-resonant 
components p^K+n^, /o(980)A'+7r^, 7r+7r~if*°, and 
7r+7r^/f+7r~; the feed-down components (1260)^^+, 
/■i:+(1270)7r-, and i^+(1400)7r-; and background compo- 
nents from qq continuum (qq), charmed _B-decays (6 
c), and charmless i?-decays (b —>■ s, u, d). 

The PDFs for the signal are separated into two cate- 
gories: correctly reconstructed events and self-cross-feed 
(SCF) events. The SCF events include at least one track 
that is taken from the accompanying B meson decay. 
For correctly reconstructed events, a sum of two Gaus- 
sians with a common mean is used for the Mbc and AE 
shapes. The M^^r and Mx-n distributions are modeled by 
relativistic Breit-Wigner functions. The p°, /2(1270) and 
K*^ resonance parameters are fixed to their PDG val- 
ues [2l|. Parameters of the /n(980) resonance shape are 
fixed to the results of Ref . [23] ; these values have higher 
precision than the corresponding PDG values. PDFs for 
the SCF components are modeled using Kernel Estima- 
tion [iil of SCF MC distributions. For the M^^ and Mk^ 
PDFs of non-resonant components, a threshold function 
and/or Chebyshev polynomials are used. The Mbc and 
AE shapes for the signal PDFs are calibrated using a 
large B^ D^tt^, K^tt~'7t~ control sample, 

to take into account small differences observed between 
MC-simulated events and data. 

The PDF shapes of the qq background are modeled 
with an ARGUS 2J] function for Afbc, linear functions 
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FIG. 2: (color online) Projection of the 4D fit results onto 
(a) Mbc, (b) AE, (c) M^rvr, (d) Mktv, with the other vari- 
ables required to satisfy (except for the variable plotted) 
Mbc € (5.27,5.29) GeV/c^ AE G (-0.045,0.045) GeV, 
M^^ e (0.62,1.04) GeV/c^ and Mk^ £ (0.84,0.94) GeV/c^ 
The curves are for the p^K'^-K~ (solid-shaded), the sum of 
p°A'*° and /o(980)iC*° (dashed), /2(1270)A'*° and the sum 
of feed-down modes (dot-dashed) , the sum of the backgrounds 
(dotted), and the total (solid). 



for AE, and combinatorial shapes for Mtt^ and Mktt- 
For b ^ c background, the PDFs are obtained separately 
for correctly reconstructed i^T*" and for random Kir com- 
binations. The fraction of each component is fixed from 
the MC simulation. The PDF shapes for b s,u,d 
background are modeled with non-parametric PDFs us- 
ing Kernel Estimation [2^. 

The following parameters are floated in the 4D fit: the 
yields of the signal modes (given in Table [l| and back- 
ground yields of & ^ c and qq; the parameters of the qq 
PDF describing the Mbc, AE and combinatorial shapes 
of Mtttt and Mk-k- The branching fractions of the feed- 
down components are fixed to the results of Ref. [25j . 
The yield of /2(1270)X*° is fixed to 43.0 events as ob- 
tained from two-dimensional A/bc-AiJ fitting in bins of 
A/tttt as discussed later. The remaining parameters are 
fixed to values obtained from MC simulations. 

The fit projections are shown in Fig. [21 and the re- 
sults are summarized in Table HI There are moderate 
correlations between some modes, which we check by 
fitting an ensemble of GEANT-simulated MC samples. 
We find a negligible effect on the measured signal yields. 
The branching fraction of each mode is determined by 
B = Y/ (eMcepiD-^Bs), where Y is the fitted signal yield, 
Emc is the event selection efficiency including daughter 
branching fractions as obtained from MC simulation and 
epiD is an efficiency correction (epiD=0.96) for PID that 
accounts for small differences between MC and data. The 
production rates of B^B^ and B'^B' pairs are assumed 
to be equal. 



TABLE I: The signal yield Y and its statistical uncertainty, 
corrected MC efficiency e (assuming /l = 0.5 for 
p^K*^), significance S including the systematic uncertainties, 
measured branching fraction B and the upper limit (UL) at 
the 90% confidence level Sul. For non-resonant decay com- 
ponents, e, B and Sul are obtained for Mktt £ (0.75, 1.20) 
GeV/c^ and M-^-,^ £ (0.55, 1.20) GeV/c^ assuming phase space 
distributions. For the branching fraction, the first (second) 
uncertainty is statistical (systematic). 



Mode 


Y 


£ 
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B 


BvL 




(events) 


(%) 


i^) 


(io-«) 


(io-«) 




77 ^+28.6 
' ' -"-27.9 


5.73 


2.7 


o 1+0.8+0.9 
^•-■--0.7-0.5 


< 3.4 


/o(980)A*o 


ri 9+20.4 


5.56 


2.5 


1 4+0.6+0.6 
-'-•*-0.5-0.4 


< 2.2 


pOA+TT- 




11.15 


5.0 


2.8 ±0.5 ±0.5 




/o(980)A+7r- 




11.43 


3.5 


1.4 ±0.4+°;^ 


< 2.1 


TT+TT^A*'^ 


200.711I-] 


6.74 


4.5 


A i: + l.l+0.9 
^•^-1.0-1.6 






_r 4+54.9 
^^•^-44. 9 


6.84 


0.0 


r, 1+1.2+1.4 
'^•J^-1.1-0.8 


< 2.1 



The fit yields the first observation for 5° — > p^K^n^ 
with a significance of 5. Oct. The significance is defined as 
■y/— 2 ln(£o / 'Cmax) , where Cq (£max) is the value of the 
likelihood function when the yield is fixed to zero (allowed 
to vary). We include systematic uncertainties by smear- 
ing the likelihood function with a Gaussian whose width 
is equal to the systematic uncertainty (discussed below). 
We also find evidence for /o(980)A'+7r^ with a 

significance of 3.5cr, and evidence for B^ 7r+7r^A'*° 
with a significance of 4.5(t. For 5° p^K*'^ and 
/o(980)Ar*", we observe excesses of events with 
significances of 2.7a and 2.5f7, respectively. For the 
non-resonant decay components, the B and smc values 
correspond to the ranges Mk-k G (0.75,1.20) GeV/c^, 
A/tttt G (0.55, 1.20) GeV/c^, and assume three-body phase 
space distributions. For modes with less than 3(t signif- 
icance, we also list a 90% confidence level (C.L.) upper 
limit. This limit is determined via 



C^C{B)dB 



= 90% . 



(2) 



/o C{B)dB 

The sources and sizes of systematic uncertainties in 



TABLE IL Summary of systematic uncertainties (%) in the 
efficiency (e) determination. 



Source 


pA* 


/A* 


pKn 


fKn 


ttttA* 


ttttAtt 


MC statistics 


±0.5 


±0.7 


±1.3 


±1.7 


±1.3 


±2.1 


Tracking 


±4.2 


±4.2 


±4.2 


±4.2 


±4.2 


±4.2 


PID 


±3.7 


±3.7 


±3.7 


±3.8 


±3.8 


±3.7 


TZqq cut 


±3.4 


±3.4 


±3.4 


±3.4 


±3.4 


±3.4 


Nbb 


±1.4 


±1.4 


±1.4 


±1.4 


±1.4 


±1.4 


./l 


+ 16.7 

-18.9 












Sum 


+ 18.0 

-20.1 


±6.7 


±6.8 


±7.0 


±6.9 


±7.0 
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TABLE III: Summary of systematic uncertainties (events) in 
the signal yield (Y) extraction. 



Source 


pK* 


fK* 


pKn 




TvnK* 


ttttA'tt 


Fitting PDFs 


+4.4 
— 5.4 


+ 12.7 
— 11.8 


+ 5.8 
— 9.1 


+ 24.1 
— 23.6 


+ 18.1 
— 17.4 


+29.4 
— 27.9 


ff 2(1270)K*'> 


+ 11.0 
— 11.3 


+5.9 
-6.4 


+0.3 
-0.3 


+0.3 
—0.1 


+ 13.9 
-13.7 


+30.0 
-35.4 


^feed-down 


+0.6 
— 1.4 


+0.1 
-0.1 


+4.7 

— 1.5 


+0.3 
—0.4 


+8.7 
— 3.8 


+3.2 
-1.9 


fb — 'S,U;d 


+ 1.9 
-2.1 


+0.1 
-0.0 


+7.0 
-9.8 


+0.3 
-0.4 


+0.0 
-1.2 


+3.7 
-0.8 


/SCF 


+2.1 
-2.1 


+ 1.2 
-1.2 


+ 19.9 
-20.6 


+ 7.4 
-7.3 


+8.2 
-8.3 


+ 11.8 
-11.4 


A'o*(1430)'' 


+29.0 
-0.0 


+ 14.7 
-0.0 


+ 16.9 
-12.4 


+0.0 
-19.3 


+0.0 
-54.8 


+69.1 
-0.0 


Fitting bias 


+2.7 
-0.0 


+4.9 
-0.0 


+ 11.2 
-0.0 


+0.0 
-10.2 


+0.0 
-26.6 


+0.0 
-29.9 


Interference 


+6.6 
-5.6 


+2.3 
-0.9 


+ 14.7 
-17.3 


+4.3 
-0.0 


+3.8 
-3.6 




Sum 


+31.5 
-12.2 


+20.5 
-12.3 


+34.8 
-31.3 


+ 25.6 
-32.9 


+35.5 
-69.8 


+76.2 

-42.6 



the efficiency determination and the yield extraction are 
summarized in Tables HIl and [TTTl respectively. The main 
sources of efficiency uncertainties are tracking (4.2%), 
PID (3.7%-3.8%), MC sample statistics (0.5%-2.1%) and 
the TZqq requirement (3.4%). Table HTl also includes the 
uncertainty from N^g (1.4%). While this does not affect 
the efficiency determination, it leads to a multiplicative 
uncertainty in B. An additional uncertainty in the effi- 
ciency for B'^ p'^K*'^ arises from the unknown fraction 
of longitudinal polarization (/l)- For our central value, 
we take /l = 0.5 and estimate the uncertainty by con- 
sidering the two extreme cases /l = and /l = 1- The 
systematic uncertainties in the yield extraction are ob- 
tained by varying all fixed parameters of the PDFs by 
±1(7, feed-down yields by ±3cr, and the fractions of SCF 
and b s,u,d backgrounds by ±50%, respectively. We 
consider the effects of higher K*'^ resonances by includ- 
ing a PDF for B° p°K^{U30f and repeating the 4D 
fit with its yield floated by extending the fitting region in 
Mxtt to 1.5 GeV/c^; the resulting changes are included 
as a systematic uncertainty. 

We study the effects of possible interference among 
/o(980), /2(1270), and non-resonant tt+tt" modes by in- 
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FIG. 3: Signal yields obtained from the two-dimensional fits 
to Mbc and AE in bins of M-mr (left) and Mktv (right) up 
to the higher-mass regions. Solid curves show the results of 
the two-dimensional binned fit, and dashed curves show the 
contributions of non-resonant ■k'^-k~ (left) and the sum of non- 
resonant K+TT- and Ko{U30f (right). The vertical lines 
show the nominal 4D fit regions. 



eluding interference terms with variable phases in the 
Mtttt relativistic Breit-Wigner function. The effect is es- 
timated by refitting with this modified PDF; the resulting 
shifts in the yields are included in the systematic uncer- 
tainties. We obtain the systematic uncertainty due to 
possible interference between K*'^, Kq{1430)'^ , and non- 
resonant K~^TT~ in the Mktt mass spectrum in a similar 
manner. Uncertainties due to possible fitting bias are de- 
termined using a large sample of MC-simulated events. 
We assign the small biases found in the MC simulation 
as systematic uncertainties. 



To verify the large contribution from non-resonant 
components (see Table U), we study background- 
subtracted Mttt^ and Mk-r spectra. These spectra are 
obtained by binning the data in M^^^^ or Mktt and, for 
each bin, fitting the two-dimensional Afbc-Ai? distribu- 
tion to determine the sum of resonant and non-resonant 
yields. Figure [3] shows these yields as a function of Mtttt 
and Mktt- Relativistic Breit-Wigner functions arc used 
as PDF's for the resonances with their parameters fixed 
to their PDG values [2l| . The PDFs for the non- resonant 
contributions are modeled by threshold functions using 
MC-simulated events. 



In summary, we have made the first observation of 
the three-body decay 5° p^K^ir^ with 5. Oct signif- 
icance and obtained the first evidence for non-resonant 
/o(980)is:+7r- and 5° 7r+7r-A'*o decays. The 
corresponding partial branching fractions are measured. 
For the 5° ^ p°K*° and 5° /o(980)A'*o modes, we 
find approximately 2.6cr signal excesses and obtain the 
results listed in Table[Tl Our result for these two-body de- 
cays are 2cr and Icr lower, respectively, than in the previ- 
ous measurement ■ We have also searched for the fully 
non-resonant four-body decay — > tt+tt^A'+tt^ and 
calculated a 90% C.L. upper limit on its partial branching 
fraction. Our results for the non-resonant modes are the 
first such results and may help us understand the polar- 
ization puzzle in pK* decays. With additional B VV 
data, these measurements can be used to constrain mod- 
els of new physics [2^ . 
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